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Abstract. We study the decay of a dark soliton in a homogeneous Bose-Einstein condensate. 
We give an analytical treatment of a decaying soliton, deriving an expression for the soliton 
velocity and decay time in the absence of thermal noise. We test the result against numerical 
simulations of a spatially confined system and find good agreement in the regime of low 
temperature (fc^T -C fi). Thermal fluctuations are found to slow the escape of the soliton, 
extending its lifetime beyond the predictions of the noise-free theory; the effect becomes 
significant at a characteristic temperature kgT ~ fi. This stabilization by noise allows us 
to infer an analytical lower bound for the dark soliton decay time. 



1. Introduction 

The observation of dark JU and bright solitons in Bose-Einstein condensates (BEC) 
has sparked much renewed interest in nonlinear excitations of one-dimensional (ID) 
superfiuids. As localized translating solutions of the repulsive, completely integrable 
nonlinear Schrodinger equation [3] 31, dark solitons also form long-lived excitations in 
trapped, finite-temperature Bose-Einstein condensates. Unlike quantum vortices |5|, dark 
solitons in simply connected geometries lack topological protection and thus the effect of 
thermal fluctuations on their stability is of particular interest. They have been created through 
a variety of experiments, including phase imprinting 0][6|-|9], using slow light ifTOll . shock 
waves ifTTl . by merging independent BECs |[T2l . or through superfiuid counterfiow [ 1 3 1 . 
Meanwhile, theoretical effort has focused on disturbing the ideal homogeneous soliton with 
thermal lfT4HT8l and quantum 1 19-25 1 fluctuations, or by breaking integrability |26| with 
a disordered potential 11271 . a stochastic perturbation 11281 . by parametrically driving the 
soliton [29 1, or confining the ID system in a non-homogeneous external potential l30rl32l . 

Significant theoretical progress has been made in the low temperature regime. In the 
zero temperature limit a Bogolyubov analysis of the breakdown of integrability due to weak 
transverse confinement predicted lifetimes very similar to experimental observations 1126) . 
Numerical simulations based on a two-fluid dynamical model of the BEC coupled to a thermal 
component lfl6l were able to reproduce the general decay phenomena observed at relatively 
low temperature JTj. More recently, application of Hamiltonian truncated Wigner theory 
revealed that fluctuations can slow the soliton below its classically predicted speed II 1 8113 31 - 
In a grand-canonical treatment of the trapped Bose gas lfl7l . soliton decay according to 
the stochastic Gross-Pitaevskii equation (SGPE) ll34l was well approximated by applying 
perturbation theory to the purely damped GPE 11351 . provided the system temperature is much 
lower than the energy scale of the soliton. This progress highlights the need for studies 
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at higher temperatures, where thermal fluctuations can significantly modify the mean-field 
dynamics of the dark soliton. 

The pervasiveness of thermal fluctuations and their importance in ID ultra-cold gas 
systems motivates our study of their role in the homogeneous dark soliton system. In 
particular, we are interested in the effect of thermal fluctuations when the system temperature 
is comparable to the soliton energy, and the departure of soliton motion from the predictions of 
perturbation theory. In the harmonically trapped system [11 711331 the soliton moves through a 
wide range of background particle densities. This complicates the analysis as the net effect of 
thermal fluctuations depends on the local background density of the condensate. By restricting 
our attention to the homogeneous system we aim to obtain a clear point of comparison for the 
effects of fluctuations. In addition to being simpler to analyze, the homogeneous regime is 
experimentally realizable in modern highly flexible experimental setups |[3"6H3"9ll . 

We apply Lagrangian perturbation theory [ 1 7. 40 -44] to find an expression for the soliton 
velocity throughout its decay in an infinite homogeneous system. This provides a simple 
quantitative basis for the intuitive picture of dark soliton decay in BECs, and a baseline for our 
grand-canonical c-field treatment based on the stochastic projected Gross-Pitaevskii equation 
(SPGPE) 843114711 . Comparing our analytical result with SPGPE simulations of a large finite 
system, we find excellent agreement at low temperature where damping dominates over noise. 
Thermal fluctuations in the SPGPE become important with increasing temperature, and their 
main effect is to slow the soliton relative to damped GPE theory, in clear departure from our 
analytical predictions. We thus find an expression for the decay time of a dark soliton that 
serves as a useful lower bound for the lifetime in the grand-canonical c-field theory. 

The structure of this article is as follows. In the remainder of this section we summarize 
the SPGPE theory for quasi-equilibrium systems, and the reduction to a ID effective theory. 
In Sec. |2] we apply Lagrangian perturbation theory to the dark soliton in a homogeneous 
background, finding analytical expressions characterizing the decaying dark soliton. In Sec. [3] 
we present numerical simulations of the dark soliton evolution, comparing our analytical 
result with the results of SPGPE and damped GPE numerical simulations. 

1.1. Stochastic projected Gross-Pitaevskii formalism 

The SPGPE provides a grand-canonical theory of finite-temperature BECs based on the 
truncated Wigner approxmation [45-49] (see also the path-integral treatment [34 50-53]). In 
addition to the usual GPE terms associated with Hamiltonian evolution, the equation of motion 
contains damping and noise terms describing interactions of low-energy atoms with a high- 
energy thermal reservoir. By incorporating the energy cutoff defining the thermal reservoir 
explicitly, the stochastic PGPE has the advantage that it may be consistently formulated in 
1, 2 and 3 spatial dimensions, is able to quantitatively model experiments, and is valid right 
across the phase transition ll54l . 

A finite temperature, partially condensed Bose gas consists of a BEC, a number of low- 
energy excitations, and a large number of high-energy modes that are weakly occupied. To 
arrive at the SPGPE, the low-energy region beneath an energy cutoff e cut is described by 
a classical field (C-field) theory ll47l . and we refer to this region as the coherent region. 
The incoherent (I) region of phase-space (above e cut ) is coupled to the coherent region via 
non-Hamiltonian two-particle scattering processes which take the form of damping and noise 
terms in the equation of motion. Neglecting the less important particle conserving processes 
(scattering terms of Ref. 0461 which are small near equilibrium), the SPGPE reduces to the 
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simple growth SPGPE |47ll49l 

<#(x, t) = V | ^(i - 7) - £gp) (x, t)<ft + dW 7 (x, t) J . (1) 

Here the GP-operator is L GP i> = (-h 2 V 2 /2m + V ext (x) + U \tp\ 2 )ip, and U = 47rft 2 a/m 
for s-wave scattering length a. The projection operator V imposes an energy cutoff, formally 
eliminating modes above e cut as part of a thermal reservoir described by \i, and T. The noise 
is a complex-gaussian Wiener-process,with non-vanishing correlation function, 

(dW;(^,t)dW-y(^',t')) = ^^S{x - x')S(t - t')dt. (2) 

This form of the noise is obtained from a first-principles microscopic treatment of reservoir 
interactions [46 49], and ensures that equilibrium ensemble properties are independent of 7. 

The dimensionless dissipation rate of the stochastic-PGPE theory can be calculated near 
equilibrium in 3D systems (49), and takes the form 



7 = 7oV ~-M— h $ 



e 



/3p(fc+l) 



P 2/3e cut k 

k=l 



l,fc 



2 



(3) 



where 70 = 4ma 2 kBT/irH 2 , and $ is the Lerch transcendent. The remarkable feature of this 
expression is that it is independent of position. It is valid in the region V ex t (x) < 2e cu t / 3 l49l , 
which typically encloses the majority of highly occupied modes, and becomes weakly 
spatially varying outside this region. The simplicity of this expression is a consequence of 
retaining the ultra-violet energy cutoff e cut in the theor>[J In general 7 is a small parameter 
depending on the system geometry, temperature, interaction strength, and particle number. 
The general phenomena considered in the present work require, as is always the case for the 
degenerate Bose gas, that 7 <C 1. 

1.2. Effective ID regime 



We consider a trap of the form 14 x t(x) = ^■LU 2 L r 2 _ + V(x) where r± = yj y 2 + z 2 . We 
identify a ID SPGPE regime defined by 

h 2 

<^lJ<huj ± < k B T. (4) 

The first condition ensures that scattering remains 3D, and the remaining conditions can be 
written in terms of length scales as 



^ < J— < J^, (5) 
mtu >]_ V rnuj± y mui']_ 

which imposes a reduction of the BEC to ID, while preserving the 3D character of the thermal 
cloud. We can thus obtain a ID reduction of the SPGPE with effective interaction strength 
g = 2hui±a, yet retain the damping rate given by the universal form (0. 

Formally eliminating the transverse degrees of freedom, we obtain the ID SPGPE 

1 

Ti 

where now L GP ijj = (-h 2 d 2 /2m + V(x) + g\^\ 2 )ip, and 

(dW*(x,t)dWJx' ,t)) = 2lk f T 8{x - x')dt. (7) 

' h 

X A consistent choice of the ultra-violet cutoff is typically of order e cu t < 3fi 1471 . 



di/j(x, t) = V { -(i - 7) Gu - L GP ) tp(x, t)dt + dW 7 (x, t) \ . (6) 
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For our analytical treatment it is convenient to work in dimensionless units; we also set 
V(x) = 0. We neglect the projector (by setting V = 1), and rescale in units of healing length 
£ 2 = h 2 /mu and speed of sound c = ^//jjrn, so that \ — x /£> T — C V£ aie dimensionless 
space and time variables. Rescaling the wavefunction ip = s/nu, with background density 
| 2 = n = u/g, we arrive at the dimensionless form of the SGPE 

du 1 d 2 u ,. t0 

with perturbation 

lR{u) = -*7 + (M 2 - + £(x,r). (9) 

The non-vanishing correlation of the Langevin noise is 

(C(X, t)SU, r)) - 2jaS( X - X ')S(r - r'). (10) 

The size of the rescaled temperature a = ksT/ p, compared to unity thus plays a central role 
in any nonlinear phenomena for finite temperature atomic superfluids. 



2. Dissipative dynamics 

To gain some understand of the different aspects of dissipation in soliton decay, we first 
consider the evolution in the absence of thermal noise. Our primary concern then is to 
understand effect of the terms explicitly proportional to 7 in Eq. © on a dark soliton that 
would otherwise propagate through a homogeneous background without change to its depth 
and velocity. 

The dynamical problem of the soliton motion may be solved by using a collective 
variable ansatz for the wavefunction. Assuming a specific form for the wavefunction 
u(x, 0,1, ■ ■ ■ , a n ), where the collective variables at are functions of time only, we minimize 
the action within the subspace of such wavefunctions. This produces a set of Euler-Lagrange 
equations for the variables aj, which may admit an analytic solution. In contrast to the 
trapped system, the homogeneous system has a subtle feature that impedes a direct approach to 
perturbation theory, namely, the non-vanishing density for £ <C |x|. We must instead consider 
a regularized Lagrangian density [40] for the homogeneous NLSE [Eq. ((HJ with jR(u) = 0], 
found by appealing to Nothers theorem: 



„ i f *du du* \ ( 1 \ 1 



du 2 



^(H 2 -i) 2 . (ID 



dx 

The (1 — l/|u| 2 ) factor regularizes the Lagrangian density to give a vanishing contribution as 
X — > ±00, and thus provides a suitable, divergence-free, starting point for a perturbative 
treatment of the dissipative dark soliton. By minimizing the action of the Lagrangian 
density ( fTTT i for the ansatz u = u(x> Q>i, ■ ■ ■ , o, n ), while allowing for a non-zero value of 
the perturbation, jR(u), one finds the Euler-Lagrange equations 

(/->•<)■ 

where L = J_ C dx, and denotes the real part. For a given perturbation and choice of 
ansatz, ( flZb provides equations of motion for the collective variables describing the system. 
The wavefunction of a dark soliton in a homogeneous system is 

^ = ^(il + l tanh (^l\] e -^/n, (13) 
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where = £/yl — v 2 /c 2 is the soliton width. Rescaling as above and dropping the rotating 
phase factor [consistent with {§)], the soliton wavefunction can be written as u(\, t) = 
iv + fttanh[ft(x — vr)) , with v = v/c, and v 2 + k 2 = 1. The simplest ansatz for the 
perturbed soliton is then 

u(x, r) = k(t) tanh [6{t)( X - x s {t))] + w(t), (14) 
with constraint v 2 + k 2 = 1. The functions ft(r), <5(r) and x s (t) describe the soliton 
depth, inverse-width, and location. The Lagrangian arising from this ansatz is L = 
2dx s /d,T [tan -1 (ft/V) — vk\ — (2/3) (k 2 <5 + k 4 /<5). Solving the Euler-Lagrange equations 
(fT2l without the effects of dissipation and noise then gives 5 = n, dn/dr = 0, and 
dx s /dT = u, as expected. 



2.1. Soliton velocity 

We now solve for the evolution of the soliton velocity with dissipation. We consider the effects 
of the perturbation (|9j while neglecting the noise (setting dw(x, t) = 0), so that it reduces to 

/ 1 Q 2 U \ 

1 R(u) = -i 1 [--— + (\u\ 2 -l)u). (15) 



2d X 2 

We can expect this term to evolve a soliton with initial speed i>i = Vi/cX0Vf = 1. Variational 
results for solitons under the effect of a number of types of dissipation have been obtained 
previously (see for example the review articles [3, 42], or [40 40 431). However, an analytic 
treatment of dark soliton decay via the damped GPE defined by ([H) and (TT3T > has not been 
previously provided. 

The equations of motion for a single soliton in our dissipative system are given by ( flZb 

as 

4x£ = ± (kS+^-), (16) 



5 3 

Ai>k = -h—, (17) 
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Figure 2. (Color online) Energy and momentum (inset) as a function of time for a decaying 
dark soliton. The initial velocity is Vi = 0.3c. 



and k 2 — S 2 . Eq. ( [ToT l implies that dx s /dr = ±v, with k = ±<5 respectively, corresponding 
to propagation in the positive or negative directions. Equation ( TP71 ), with v = —nh/v, gives 
the soliton equation of motion 



d{v 2 ) 4 7 



dr 3 

with fixed points at v = and v = ±1. Given initial speed < i>$ < 1, the solution is 



v{t) = 



^(l-„2)e-^ + . 



(18) 



(19) 



Clearly v —> vi as r — >• 0, and i/ — > 1 as r — > oo, with characteristic timescale ~ 7 _1 > 
or, in dimensional units, ~ = r<j. This result can be used to give an analytical 

description of the dark soliton evolution due to the effects of thermal damping, including the 
decay time and a number of other useful properties, providing a point of comparison for full 
stochastic simulations. As shown in Figure [TJ a dark soliton subject to damping accelerates 
from 0<fi<ltoi/ = l, corresponding to a particle with negative effective mass. 

An immediate consequence of obtaining a damping solution is that we can easily evaluate 
the dark soliton energy E = (4/i/3)(l — v 2 ) 3 ! 2 , mass m = — (4/i/c 2 )v / l — v 2 , and 



momentum Ic/fi = —2v\[\ — 



E(r) 
m(r) 



Ei 



2 tan" 1 (Vl -v 2 /v): 

l + ^(e 4 ^/ 3 -l) 
e 3~ 



4~fT 

3 



J(r)c _ -2u iy /T 



-2 7 r/3 



(l_ |y 2) e -4 T r/3 + i 



2 tan" 



-2 7 t/3 



(20) 



(21) 



(22) 



An example of the soliton mass is shown in Figure Q] (inset), and the energy and momentum 
are shown in Figure |2] 
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Figure 3. (Color online) The decay time [Eq. 1231 1 of a damped soliton in a homogenous 
BEC, versus initial velocity Vi, for a soliton that is resolvable up to a velocity vt = 0.9c. 



We can also define a decay-time Td for the soliton as the time for the velocity to evolve 
between v± and vt. Using ( fT9] l we find 



The decay-time diverges for Vi — > or Vf — > 1, but the two divergences have a different 
physical origin. When v t — > the soliton becomes immune to damping. The state v = is an 
eigenstate of parity. In practice any small thermal or quantum fluctuation will cause symmetry 
breaking and destabilize the v- t = state, forcing a direction of motion to be selected. The 
divergence as Vf — > 1 is due to the asymptotic nature of the decay. As the soliton evolves it 
will eventually become too shallow to be experimentally resolved, at which point the soliton 
has decayed for all practical purposes. The dependence of soliton decay time on z/j is shown in 
Figure[3]for a range of damping rates and v$ = 0.9. The logarithmic dependence on velocities 
seen in (l23l is qualitatively similar to that found by Fedichev et al (Eq. (13) of Ref. ifPfll ). via 
different reasoning. 

3. Stochastic simulations 

Our primary interest in this work is the influence of noise on the soliton velocity. In 
this section we numerically simulate the evolution of a dark soliton in a large, finite, one 
dimensional system, in order to assess the effect of thermal fluctuations on the soliton's 
velocity during its evolution. Our dimensionless stochastic equation of motion is 



with Lgpu = [—(l/2)d 2 /dx 2 + \u\ 2 + V(x)]u, and complex-Gaussian noise defined by 
non-vanishing correlator 




(23) 




(24) 



(dw*( X , T)dw( X ',T')) = 2 1 a5(x' - x)*(t - r')dr. 



(25) 
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Figure 4. A single soliton imprinted on a sample from the SGPE ensemble. The external 
potential is zero within the condensate, and provides steep confinement near \ = ±60, The 
temperature is a = 0.3. 
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Figure 5. (Color online) Ensemble-averaged velocity (extracted from soliton depth) of a dark 
soliton according to the SPGPE 124 \ at temperature is a = 0.1. Also shown are the results 
from the damped GPE, and the analytical solution i\9\ . The horizontal line shows the highest 
numerically resolvable velocity (which differs from 1 due to thermal fluctuations). 



The potential is chosen as V(x) = for —60 < \ < 60, with steep confinement at the 
endpoints, creating a homogeneous region of length Lh ~ 120£. We begin by creating 
a finite temperature equilibrium ensemble of trajectories for the confined system using the 
SPGPE. We initialize our dynamical simulations by imprinting a positively directed soliton 
wavefunction of given depth onto each trajectory. We use = 0.3c, and interaction parameter 
g = 0.02 for all of our SPGPE simulations. For reasons of numerical expedience we also use 
7 = 0.05. A detailed treatment of specific experimental systems requires that we include 
the full temperature dependence of 7. In this work we focus on the role of fluctuations, for 
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Figure 6. (Color online) Velocity of a dark soliton according to the SPGPE I24t at temperature 
cr = 0.3. All other parameters are the same as in Fig. [5] 



fixed damping^. From Figure[3]we can expect the soliton to damp out over a distance of order 
~ 40£. We choose the initial position of the soliton as x s (0) = —30 so that it remains well 
inside the homogeneous region throughout. Figure|4]shows an element of the initial ensemble 
after the soliton has been imprinted. 

We simulate the decay of a single soliton for 400 SGPE trajectories, and compare the 
results with damped PGPE simulations of a soliton with the same initial speed. To calculate 
the velocity of the soliton numerically, the minimum density of the condensate is found within 
the region |x| < 3Lh/8, and the corresponding soliton velocity is calculated. Eventually the 
soliton decays sufficiently that it becomes indistinguishable from the background fluctuations. 
After this point the soliton is unresolvable, so we do not show soliton velocities after this time. 

In figures[5][6] and|7]we compare the soliton velocity from the analytical result ( fT9l , with 
numerical simulations of the damped PGPE [© and (fT3Tll, and with ensemble averages from 
d24l . The comparison is made for temperatures a = 0.1, 0.3 and 0.6. At low temperature 
[Figure [5], we see excellent agreement between the average SPGPE result and perturbation 
theory. At the higher temperatures the SPGPE shows a significant departure from both 
damped PGPE and perturbation theory. As temperature increases there is also a significant 
increase in the variance of the soliton velocity. Most notably, the stochastic fluctuations have 
the principal effect of slowing the decay process: the average soliton velocity in SPGPE is 
always lower than those found from damped GPE or soliton perturbation theory. 

Recent work on soliton core-filling [ 18 1 indicates that phase fluctuations in an ensemble 
of SPGPE trajectories should lower the expectation value of the soliton density and hence 
its instantaneous ensemble-averaged velocity. This expectation is confirmed by our analysis, 
and is a consequence of the classical relationship between the phase jump across the soliton 
(4> c ) and its central density (n s ). In a system with background density rib, the equation is 
n s = rib cos 2 (<p c /2). Thermal fluctuations generate a temperature dependent variance of 
4> c around the mean field solution. A higher temperature will thus yield a smaller value for 
(n s ) oc (cos 2 (0 c /2)). 

Based on our observations of the stochastic slowing down of the a dark soliton, we infer 



§ However note that this choice of 7 is of the order of that realized in recent experiments 11471 . 
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Figure 7. (Color online) Velocity of a dark soliton according to the SPGPE I24t at temperature 
a = 0.6. All other parameters are the same as in Fig. [5] 

that ([23} provides an analytical lower bound for the mean decay time of the dark soliton. 
4. Conclusions 

In this work we have studied a novel regime for dark soliton dynamics in finite-temperature 
superfiuids; this regime is characterized by temperatures comparable to the soliton energy, 
kgT ~ fj„ We have established that in this regime the dark soliton exhibits an interesting 
physical phenomena: noise-induced slowing-down. 

For the specific damping term arising in the high temperature classical field theory (the 
stochastic projected Gross-Pitaevskii equation, or SPGPE), we obtained an explicit solution 
for the velocity, and hence for the energy, mass and momentum, and decay time of a decaying 
dark soliton in a homogeneous BEC. Parity symmetry of the black soliton restricts these 
results to grey soliton initial conditions. We demonstrated the accuracy of the results at low 
temperature (fc^T -c /it) by comparing the soliton velocity with numerical simulations of the 
SPGPE, finding excellent agreement between the two approaches: in this regime the damping 
term dominates over the noise. At high temperature (fc^T ~ /i), the noise significantly 
modifies the soliton dynamics, and the analytical expression (T% departs from numerical 
simulations of (|24| |. However, since the noise consistently slows the decay of the soliton, our 
expression for the lifetime based on damping [Eq. (TfJt l provides a lower bound for the dark 
soliton decay time. 

We conclude that a description of dissipation that includes damping and the associated 
noise, leads to a dark soliton with increased stability. This contrasts strongly with the reduced 
stability of a quantized vortex in a spherically trapped 3D Bose gas ll55l . but is consistent with 
a recent low-temperature treatment of the dark soliton [fT8l . 

Future work will focus on more detailed modeling of specific experimental systems 
and providing an analytical treatment of the noise. The latter requires treating noise in the 
collective variable equations of Lagrangian perturbation theory — evidently a challenging task 
for the dark soliton. 
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